From percolating to dense random stick networks: Conductivity model investigation 
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In a Monte Carlo study the conductivity of two-dimensional random stick systems is investigated from the 
percolation threshold up to ten times the percolation threshold density. We propose a model explicitly depend- 
ing on the stick density and junction-to-stick conductance ratio. The model describes the transition from the 
conductivity determined by the structure of a percolating cluster to the conductivity of the dense random stick 
networks. The model is motivated by the observed densities of the sticks and contacts involved in the current 
flow. The finite-size scaling effects are also included in the description. The derived model for conductivity 
should be broadly applicable to the random networks of the rodlike particles. 
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I. INTRODUCTION 

Recently there has been an increasing interest in the net- 
works of randomly distributed stick (rodlike) particles «, due 
to the development in the area of conductive nanoparticles, 
such as the carbon nanotubes, silicon, copper, and silver 
nanowires 5 , and the promising applications in electronics 6 , 
optoelectronics 7 , and sensors 8 . The conventional electronic 
composites containing stick particles as a filler can be used 
as a conducting channel in the thin-film network configura- 
tions^—. Probably, the most important characteristics of the 
stick networks, are the conductivity, either electric 12 or ther- 
mal^. The conductivity dependence on the stick density and 
system geometry needs to be taken into account in any de- 
vice design 14 . The percolation models 1516 are often used to 
model an onset of the high electrical conductivity in the com- 
posites consisting of the conductive sticks in the insulating 
matrice a 1,7,17 . 

The percolation theory predicts that the electrical conduc- 
tivity of the composite materials with the conductive filler 
density n above, but close to the percolation threshold n c , in- 
creases with the density by a power scaling law a ~ (n— n c )*, 
with the universal conductivity exponent t ~ 1.29 for two- 
dimensional (2D) systems 15 . While the conductivity scaling 
law is expected to be applicable only near the percolation 
threshold, in many experiments the scaling law was used over 
a much larger range of concentration, but with the nonuni- 
versal values of the conductivity exponen t 7 ' 17 ' 18 . Hu et al.— 
obtained the nonuniversal value 1.5 for the conductivity ex- 
ponent using the conductivity scaling law for fitting the ex- 
perimental data for ultrathin carbon nanotube networks oper- 
ating from the percolation threshold up to about ten times the 
percolation threshold density. Several numerical studies con- 
firmed the observed nonuniversality of the conductivity ex- 
ponent when the stick density is well above the percolation 
threshold^—. Keblinski et al. 19 demonstrated that the uni- 
versal power law holds from the percolation threshold n c , to 
about twice its value 2n c . For higher stick density, n > 2n c , 
and in two limiting cases they observed that the conductiv- 
ity scaling exponent becomes: (i) slightly higher than 1 when 
the junctions are superconductive and only the stick conduc- 
tance is the limiting factor for the current flow through the 
system and (ii) close to 1.75 when the sticks are supercon- 



ductive and the contact conductance is the limiting factor. Li 
et al. 20 showed that the conductivity exponent significantly 
varies with the junction-to-stick conductance ratio for lower 
stick densities up to 2n c . The broad range applicability of the 
conductivity scaling law was explained by the presence of the 
long-range correlations in the distribution of the conductive 
sticks in the system 22 . We will demonstrate that the nonuni- 
versality of the conductivity exponents is a consequence of a 
transition from the percolating to dense stick networks. 

In this paper, we numerically investigate the conductivity 
of the stick systems from the percolation threshold up to ten 
times the percolation threshold density. We show that it is not 
appropriate to use a simple scaling law to describe the conduc- 
tivity dependence on the density both for finite and dense sys- 
tems. Based on the Monte Carlo simulation results, a model 
is proposed describing the conductivity dependence on the 
stick density and the different junction-to-stick conductance 
ratios. The proposed model is valid for the different stick- like 
nanoparticles (e.g., the carbon nanotubes and the nanowires). 
The model is motivated by the observed structural character- 
istic (i.e., the density of the total sticks and contacts involved 
in the current flow through the system). The finite-size ef- 
fects, especially pronounced in the vicinity of the percolation 
threshold, are included in the generic description for the con- 
ductivity of stick systems. 



II. NUMERICAL METHOD 

Monte Carlo (MC) simulations are coupled with an efficient 
iterative algorithm implemented on the grid platform and used 
to investigate the conductivity of stick systems 2Q ' 23 i 24 . We 
have considered the two-dimensional systems with isotropi- 
cally placed widthless sticks of length I. The centers of the 
sticks are randomly positioned and oriented inside the square 
system with size L. Two electrodes (i.e., conducting bars) 
are placed at the left and right sides. The top and bottom 
boundaries of the system are free and nonconducting. The 
free boundary conditions are more consistent with the finite - 
size rodlike nanoparticle networks in practice 711 . Two sticks 
lie in the same cluster if they intersect. The system perco- 
lates (conduct) if the electrodes are connected with the same 
cluster. The behavior of the stick percolation is studied in 
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terms of the stick density n = N/ (L/l) 2 , where N is the to- 
tal number of sticks and L/l is the normalized system size. 
The percolation threshold of the infinite system is defined by 
the critical density n c rs 5.63726 (Refsi 2 ^ 6 .). To evaluate 
the conductivity of the stick systems we introduce two dif- 
ferent conductances: (1) the conductance of the entire stick 
G s and (2) the conductance due to the stick-to-stick junction 
Gj. We assume diffusive electrical transport through the stick 
typical for the rodlike nanostructures (carbon nanotubes and 
nanowires) whose length is larger than the mean free path of 
the electrons 27 . According to the diffusive electrical transport 
the electrical resistance of a stick segment is proportional to 
the length of the segment 28 . In our simulations, each stick- 
stick junction is modeled by an effective contact conductance 
regardless of the type of the junction, following the simpli- 
fied approach of the authors of Refs i 6,14 ' 20 . Therefore, if two 
sticks intersect a junction with the fixed conductance Gj is 
created at the intersection point. If a stick intersects an elec- 
trode the potential of the electrode is applied to the intersec- 
tion point. Kirchhoff's current law was used to balance the 
current flow through each node of the created network. An 
iterative equation solver (i.e., conjugate gradient method with 
Jacobi preconditioner) has been employed to solve a large sys- 
tem of the linear equations following from the current la w 20 ' 29 . 
After obtaining the total current / under an applied voltage V 
the macroscopic electrical conductivity of the system is eval- 
uated as a = I/V (Ref. 30 ). Monte Carlo simulations have 
been performed for a wide range (i.e., Gj/G s = 0.001 to 
1000) of junction-to-stick conductance ratios (cf. Refs<2i~— ). 
Finally, for each set of the system parameters, the electrical 
conductivity is averaged over the Nmc independent MC re- 
alizations. To obtain the same precision for the finite-size 
systems Nmc — 64000 realizations are used for the systems 
with normalized size L/l = 10 down to Nmc = 4000 for the 
largest system L/l = 40 studied. Using the appropriate func- 
tions for the fitting data and the least-squares fitting method- 
ology 26 , good fits with high correlation factors (R 2 > 0.998) 
were obtained for all analyzed systems. 

III. RESULTS AND DISCUSSION 

As already mentioned, the numerical estimates of the con- 
ductivity exponent t are based on the linear fit of the MC 
results for the logarithms of the conductivity a and density 
n — n c (Refsii^r— ). The estimates therefore rely on the as- 
sumption that a obeys the simple power-law dependence over 
a quite extended density range. As there exists no justifica- 
tion of such an assumption, we have investigated in detail the 
sbehavior of the conductivity a as we move away from the 
critical point. A local (density dependent) conductivity expo- 
nent is defined as t(n) by^l 



The dependence of the local conductivity exponent t (n) on the 
stick density n and the ratio of the stick-stick junction conduc- 
tance (Gj) to stick conductance (G s ) (i.e., Gj/G s ) is shown 



in Fig. Q] As one can see from a coarse observation, when 
the stick density approaches the percolation threshold n c from 
above the local conductivity exponent converges to the univer- 
sal value for 2D systems t(n c ) ps 1.29 for all Gj/G s values. 
The fine behavior of the local conductivity exponent for finite- 
size systems in the vicinity of the percolation threshold will be 
discussed later in this section. With the increasing concentra- 
tion n, the local conductivity exponents t(n) change quickly 
from the universal value t(n c ), taking the values in a wide 
range 1 < t(n) < 2. From Fig. [1] one can see that the local 
conductivity exponents t(n) for the conductance ratio higher 
than 2 (Gj/G s > 2) is a monotonically decreasing function of 
the stick density n which converges to 1 from above. Some- 
what surprisingly, for the conductance ratios lower than 1 (i.e., 
Gj/Gs < 1), the local exponent t(n) is not a monotonic func- 
tion and has a local maximum. The observed density where 
the local conductivity exponent reaches a maximum is de- 
creasing with the conductance ratio G , / G.. 
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FIG. 1 . (Color online) The dependence of the local conductivity ex- 
ponent t{p) on the stick density n and junction-to-stick conductance 
ratio Gj /G s - The points are MC simulation results obtained using 
Eq. Q} for the system size L/l — 20. The values are given for the 
conductance ratios G 3 /G s = 0.001, 0.01, 0.1, 0.2, 0.5, 1 (filled), 
and their inverse values 1000, 100, 10, 5, 2 (transparent). The error 
bars are smaller than the size of the points. The star marker denotes 
the expected universal value for the conductivity exponent at the per- 
colation threshold i(n c ). The lines represent the local conductivity 
exponents t(n) obtained from the conductivity model for an infinite- 
size system, Eq. (3)- 

To explain the observed behavior of the exponent t(n) at 
the higher densities n > 2n c , one needs to look into the struc- 
ture of the dense conducting stick systems. Figure|2]shows the 
densities of the sticks n 1 and junctions that carry the cur- 
rent through the system. For sufficiently high stick densities 
(n > 2n c ), one can see that almost all the sticks and junctions 
in the system contribute to the conductivity and that the den- 
sity of the current-carrying junctions increases with the stick 
density n by a square power law ~ n 2 . The reason for this 
is that the mean number of contacts per stick is proportional 
to the stick density, see Ref. 39 . Also, for a sufficiently high 
stick density n the current-carrying stick density n 1 is propor- 
tional to n. Therefore, when the stick density is well above 
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the percolation threshold n c , the conductivity of the system 
can be modeled as an equivalent serial conductance n sticks 
in parallel and n 2 junctions in parallel 



1 



bn^/Gs+n^/Gj 1 



(2) 



where b is a constant parameter. One can see that the square 
term n~ 2 /Gj, originating from the junctions, converges faster 
to zero than the linear term bnr 1 /G s . This explains the con- 
ductivity exponent t(n) approaching to 1 when the stick den- 
sity is sufficiently high (i.e., n ^> G s /Gj) and the existence of 
the local exponent maximum in Fig.Q] If the sticks are much 
more conductive than the junctions (e.g., Gj/G s = 0.01) the 
density where the local conductivity exponent starts to con- 
verge to 1 is high and computationally unreachable in the MC 
simulations shown in Fig. Q] Only in the limiting case when 
the sticks are superconductive and the conductance ratio ap- 
proaches zero, i.e., Gj/G s — > 0, should the conductivity ex- 
ponent t(n) converge to 2 with the increasing density n, which 
is consistent with Keblinski et alJ^. In the other limit, when 
the junctions are superconductive (i.e., Gj/G s —> oo) the lo- 
cal conductivity exponent t(n) should have the fastest conver- 
gence to 1. 
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FIG. 2. (Color online) The density of junctions nj and sticks n 1 
involved in the current flow through the system is compared with 
the density of all junctions rij and sticks n in the system of size 
L/l — 20. For higher stick densities n almost all junctions and sticks 
will carry some current. The error bars are smaller than the size of 
the points. Inset: The density ratio of the current-carrying junctions 
to current-carrying sticks nj /n 1 is higher than the density ratio of 
all junctions to all sticks n,/n in the system. At the percolation 
threshold this ratio is about 2 [i.e., n^/n = 2.0(1)]. 

At the densities close to the percolation threshold n c , only a 
fraction but not all the sticks and junctions in the system con- 
tribute to the conductivity, by carrying some current. From 
Fig. |2] (inset), one can see that at the percolation threshold 
n c , the density of the current-carrying junctions is about two 
times higher than the density of the current-carrying sticks 
[i.e, rij/n 1 — 2.0(1)]. From the framework of the perco- 
lation theory we cannot determine a density-dependent fac- 
tor of proportionality in the conductivity power law [i.e., 



(7 ~ (n — n c )*]. Instead, we fit the factor of proportion- 
ality with an expression for the dense systems [i.e., Eq. (f2)] 
and obtain 1/ [bn^/Gg + (n + n^ 2 /Gj] . This relation 
explicitly includes the previous observation that there is al- 
most exactly two times more current-carrying junctions than 
current-carrying sticks at the percolation threshold. For a gen- 
eral conductivity description of the infinite-size systems we 
obtain 



a = a 



(n - n c y 



&rc*-7G s + (n + n c )*- 2 /<V 



(3) 



where a — 0.027(1) and b = 0.061(3) are fitting parameters 
calculated using the least-squares fitting methods. The solid 
lines in Fig. Q] denote the local conductivity exponents t(n) 
calculated from Eq. (fl}, using the model for an infinite sys- 
tem given by Eq. OJ, for a wide range of conductance ratios 
Gj/G s = 0.001 to 1000. Deviations between the modeled 
and MC values for local conductivity exponent t(n) are com- 
parable to the statistical errors. 




FIG. 3. (Color online) Simulated current (a) and (b) at different sys- 
tem sizes and (c) and (d) junction-to-stick conductance ratios. The 
current through a stick / is given relative to the maximal current in 
the system / max . There is a large difference in the fraction of the sys- 
tem involved in the current flow between the two nominally identical 
films in term of density (n = n c ) and junction-to-stick conductance 
ratio (Gj/Gs = 1) for two different system sizes L/l — 10 and 
40. The current redistribution with the increasing junction-to-stick 
conductance ratio Gj / G s is visible from (c) and (d). (c) If junctions 
are weakly conductive (i.e., Gj/G 3 — 0.01) the maximal current is 
flows along the shortest path with the least junctions, (d) For high 
junction conductance values (i.e. Gj/G s = 100), the total current 
is evenly carried by the larger number of shortest paths connecting 
electrodes. This effect is only visible at higher densities (e.g., n = 8) 
where several paths connecting electrodes exist. 
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FIG. 4. (Color online) The local conductivity exponents t(n) for 
the stick systems with increasing size L/l — 10, 20, and 40 and for 
three conductance ratio values (a) Gj/G s =0.01, (b) 1, and (c) 100. 
The direction of the increase of L/l is indicated on the graphs. The 
points are obtained from the MC simulations and calculated using 
Eq. |QJ. The error bars are smaller than the size of the points. The 
solid line represents the local conductivity exponent t(n) for the in- 
finite system obtained from Eq. ((3}, while the dashed lines denote 
the local conductivity exponents t(n) obtained from the model that 
includes finite-size effects, Eq. {7J. The star marker denotes the ex- 
pected value for the conductivity exponent of the infinite system at 
the percolation threshold n c . 



Figure [3] illustrates the structure of the percolating cluster 
[Figs. [3 a) and Ob)] and the redistribution of the current in 
the dense stick networks due to the junction-to-stick conduc- 
tance ratio increase [Figs.|3jc) and[3jd)]. The current through 
a stick / is given relative to the maximal current in the system 
Imax- As one can see from Figs. [2 a) and[3jb), the percolating 
cluster consists of a few sub-clusters connected by high cur- 
rent links. This explains why on average more junctions than 
sticks are needed to shortcut the electrodes. For a large, but 
finite-size systems at the percolation threshold, the density of 
the current-carrying junctions decreases as ~ ( y L/l)~^/ 1 ' 
with normalized system size L/l, where f3 = 5/36 for 2D 
systems 15 . Also, the density of the current-carrying sticks at 
the percolation threshold is n ~ (L/l)-P/ v . As a result 
the densities of the current-carrying sticks and junctions de- 
crease with system size [cf. Figs. [3l a ) (L/l — 10) and[3|b) 
(L/l = 40)]. Furthermore, the density ratio rij/n 1 at the per- 
colation threshold converges to a constant value with the in- 
crease of the system size, see Fig. [2] At higher stick densities 
(i.e., n = 8) one can see that current flows along many parallel 
paths connecting electrodes. An increase of the junction-to- 
stick conductance ratio Gj / G s results in the more uniform re- 
distribution of the current [cf. Figs. |3]c) (Gj/G s = 0.01) and 
[3jd) (Gj/Gs — 100)]. For weakly conductive junctions (i.e., 
low conductance ratio Gj/G s = 0.01), most of the current 
flows through a shortest path with the least junctions along. 
With the increase of the junction conductance several parallel 
paths become visible. As a result, the total current through 
the system is more evenly distributed, resulting in the higher 
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FIG. 5. (Color online) (a) Conductivity as a function of (n — n c )/n c 
is obtained from the MC simulations (points) for the stick system 
of size L/l = 20 and the junction-to-stick conductance ratio from 
Gj /G a = 0.001 to 1000 (from bottom to top). The lines denote val- 
ues obtained from the conductivity model for the finite-size systems 
given by Eq. Q. (b) The conductivity ratio between the MC simula- 
tion results ctmc and the values obtained from the model (Tmodei for 
corresponding finite-size systems, Eq. ((T). The error bars are smaller 
than the size of the points. 



conductivity. This is also expected from Eq. (0. 

If we compare the infinite system model prediction and MC 
simulation results in Fig.[T]close to the percolation threshold, 
we observe a deviation between the predicted and simulated 
values. This deviation is a result of the finite-size effects, since 
the MC results in Fig.Q]are calculated for the large but finite- 
size system (i.e, L/l — 20). The convergence of the local con- 
ductivity exponents with the increasing system size is shown 
in Fig. [4] The points are MC simulation results for the sys- 
tems with sizes L/l — 10, 20, and 40 and the solid line de- 
notes the model for an infinite system given by Eq. (01. For 
the finite-size systems close to the percolation threshold we 
observe a large deviation of the local conductivity exponent 
t(n) from the model. The local conductivity exponent de- 
creases with the decreasing system size and can be even lower 
then 1 [t(n) < 1]. This is result of a nonzero conductivity 
value for the finite-size systems at the percolation threshold 
— . Therefore, the model should be adapted for the finite-size 
systems. The finite-size scaling argument s 15 i 26 i 40 suggest that 
the conductivity a depends on the system size L as 



(n-n c )*/ 



L 



(4) 



where £(n) ~ l\n — n c \~ v is the correlation length that mea- 
sures the linear extent of the largest cluster. For 2D sys- 
tems the correlation-length exponent is v = 4/3 (RefJ£). 
For the infinite system above the percolation threshold [i.e., 
£(n)/L — > 0] the conductivity follows the simple scaling law 
and finite-size scaling function / [£ (n) / L] converges to a con- 
stant value. In the other limit, for the finite-size systems at 
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the percolation threshold [i.e. £(n)/L — > oo], conductivity 
has a nonzero value g ~ {L/l)- 1 '" (Ref^). Therefore, the 
finite-size scaling function should have a form / [£(n)/L] ~ 
[£(n) / £] , to cancel the conductivity dependence on density 
in Eq. ©. Since the finite-size scaling function / [£(n)/L] 
above the percolation threshold is a continuous and smooth 
function 15 , we approximate it by a combination of its two lim- 
iting behaviors 



L 



''(L/l) 



-t/u 



(5) 



where c is the finite-size parameter. Inserting Eq. (|5) into 
Eq. (01 the first-order approximation of the finite-size scaling 
law for conductivity becomes 



(n-n c )* +c{L/iy t / u . 



(6) 



Finally, incorporating the finite-size effects given by Eq. (0 
into the conductivity model for an infinite-size system, 
Eq. (0, we obtain the finite-size model for conductivity 



(n-n c ) t + c{L/l)- t /'' 
'bnt-yGs + in + n.y-yGj 



(7) 



close to the percolation threshold. Finally, the MC conductiv- 
ity values normalized with the stick conductance G s and fitted 
by Eq. (0 for the systems of size L/l = 20 and conductance 
ratios from Gj/G s = 0.001 to 1000 are shown in Fig. 0a). 
For all studied values of the conductance ratio Gj /G s the con- 
ductivity obtained from the model agrees with the MC results 
over the whole range of the stick density n, see Fig. a). The 
agreement between the MC results and the model is good for 
higher stick densities (n > 2n c ) (i.e., further away from the 
percolation threshold), but not so good in the vicinity of the 
percolation threshold [cf. Fig. 0b)]. Hence, in the vicinity 
of the percolation threshold the conductivity ratio between the 
MC simulation results and the values obtained from the model 
is shown in Fig.0for different system sizes L/l = 10, 20, and 
40. For all three system sizes in Fig. 0the curves look quali- 
tatively similar. Only the density where the dense-system be- 
havior becomes dominant decreases with the system size L/l. 
To improve the agreement between the MC results and the 
model close to the percolation threshold one could consider a 
further refinement of the model to include higher-order cor- 
rection for the finite-size effect. Finally, the proposed model 
for conductivity gives a good estimate of the local conductiv- 
ity exponents, as one can see in Figs.0and0 
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FIG. 6. (Color online) The conductivity ratio between the MC sim- 
ulation results <7MC and the values obtained from the model <r mo d c i 
given by Eq. (0 for different system sizes (a) L/l — 10, (b) 20, and 
(c) 40 and for three conductance ratio values Gj/G s = 0.01, 1, and 
100. The error bars are smaller than the size of the points. 

The finite-size parameter for 2D stick systems c = 2.5(1) is 
calculated using the least squares fitting methods. A compar- 
ison between the MC results and the values obtained from the 
model given by Eq. (0 is shown in Fig. [4] The dashed lines in 
Fig.0denote the local conductivity exponents f (n) calculated 
from the model including finite-size effects, Eq. 0, for sys- 
tems with increasing size L/l = 10, 20, and 40 and for three 
conductance ratio values Gj/G s — 0.01, 1, and 100. We 
see that the introduction of finite-size effects in the model sig- 
nificantly improves the quantitative description of the system 



IV. CONCLUSIONS 

In this paper, we present the results of the numerical Monte 
Carlo study of the conductivity of random stick systems for 
the wide range of densities and junction-to-stick conductance 
ratios. We observe the transition from the conductivity of 
the percolating cluster to the conductivity of the dense ran- 
dom stick networks with increasing density. Three limiting 
cases are identified for the conductivity of whole system: one 
in the vicinity of the percolation threshold, and two for high 
densities when either the junctions or sticks are superconduc- 
tive. Each of these cases has a different exponent governing 
the power-law dependence of the conductivity from density 
(i.e., 1.29, 1, and 2, respectively). As result, the exponent 
can take values anywhere in the range (1,2) depending on the 
junction-to-stick conductance ratio. For finite-size systems 
the density-dependent exponent can even take values lower 
than 1 . Therefore, it is not appropriate to use a simple scaling 
law to describe the conductivity dependence on the density 
both for finite and dense systems. We instead propose a com- 
prehensive conductivity model, derived from the behavior of 
the limiting cases. We find that the proposed description gives 
a satisfactory estimation of the conductivity and the local con- 
ductivity exponent (which is related to the first derivative of 
the conductivity) over the whole range of the stick density 
values. Finite-size effects, important for many practical real- 
izations of the random conducting networks, are also included 
in the conductivity model. The presented methodology could 
be used to describe the properties of other conducting systems 
(i.e., disks, spheres, and fibers). 
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For a rectangular system of size L x x L y , where L x is the length 
of the electrodes and L y is the distance between them, the relation 
between the system conductance G and conductivity o is accord- 
ing to Ohm's law given by G — oL x /L y . In this paper the system 
is the square-shaped L x — L y — L, which implies G = a. Al- 
though in the square system conductance is equal to the system 
conductivity, for clarity we have assumed terminology in which 
the conductance is used for denoting a single conductive element 
(stick or junction) while the conductivity is related to the entire 
system. 
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